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Abstract 

We calculate the minimal energy path over the sphaleron barrier in the pre- 
sence of fermions, assuming that the fermions of a doublet are degenerate in 
mass. This allows for spherically symmetric ansatze for the fields, when the 
mixing angle dependence is neglected. While light fermions have little influence 
on the barrier, the presence of heavy fermions (Mp ~ TeV) strongly deforms the 
barrier, giving rise to additional sphalerons for very heavy fermions (Mp ~ 10 
TeV). Heavy fermions form non-topological solitons in the vacuum sector. 



Utrecht-Preprint THU-93/17 



1 



1 Introduction 



In 1976 't Hooft [1] observed that the standard model does not absolutely conserve 
baryon and lepton number due to the Adler-Bell-Jackiw anomaly. The process 't Hooft 
considered was spontaneous fermion number violation due to instanton induced tran- 
sitions. Fermion number violating tunnelling transitions between topologically distinct 
vacua might indeed be observable at high energies at future accelerators [2,3]. 

Manton considered the possibility of fermion number violation in the standard model 
from another point of view [4]. Investigating the topological structure of the configura- 
tion space of the Weinberg-Salam theory, Manton showed that there are noncontractible 
loops in configuration space, and predicted the existence of a static, unstable solution of 
the field equations, a sphaleron [5], representing the top of the energy barrier between 
topologically distinct vacua. 

At finite temperature this energy barrier between topologically distinct vacua can 
be overcome due to thermal fluctuations of the fields, and fermion number violating 
vacuum to vacuum transitions involving changes of baryon and lepton number can 
occur. The rate for such baryon number violating processes is largely determined by 
a Boltzmann factor, containing the height of the barrier at a given temperature and 
thus the energy of the sphaleron. Baryon number violation in the standard model due 
to such transitions over the barrier may be relevant for the generation of the baryon 
asymmetry of the universe [6-10]. 

While the barrier between topologically distinct vacua is traversed, the Chern- 
Simons number changes continuously from Acs in one vacuum sector to Acs + 1 in the 
neighbouring vacuum sector. At the same time one occupied fermion level crosses from 
the positive continuum to the negative continuum, leading to the change in fermion 
number. When considered in the background field approximation this level crossing 
phenomenon predicts the existence of a fermion zero mode precisely at the top of the 
barrier, at the sphaleron. For massless fermions this zero mode is known analytically 
[11-13]- 

Considering the minimal energy path over the barrier [14,15], the fermionic level 
crossing was demonstrated recently in the background field approximation for the 
Weinberg-Salam theory [16], assuming that the fermions of a doublet are degener- 
ate in mass. This assumption, violated in the standard model, allows for spherically 
symmetric ansatze for all of the fields, when the mixing angle dependence is neglected 
(which is an excellent approximation [17,18]). 

With these assumptions we here calculate the minimal energy path over the barrier 
self-consistently, varying both the fermion mass and the Higgs mass. The presence of 
the fermions influences the energy and the shape of the barrier, which need no longer 
be symmetric with respect to the sphaleron. Neither must the fermion zero mode occur 
precisely at the top of the barrier. 
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While our main concern is the study of the barrier for light fermions up to about the 
top quark mass, we also consider heavy fermions with masses in the TeV region, because 
the possibility of rapid anomalous decay was discussed recently for heavy fermions, and 
the existence of a new type of soliton was conjectured [19], which is distinct from the 
known weak non-topological soliton [19-21] representing a weak chiral soliton. Even 
for heavy fermions we employ the valence fermion approximation, neglecting all boson 
loops and the effects of the Dirac sea (as in [19]). 

We briefly review in section 2 the Weinberg-Salam lagrangian and the anomalous 
currents for vanishing mixing angle and for degenerate fermion doublets. In section 
3 we present our ansatz, obtain the energy functional and the Chern-Simons number 
and derive the equations of motion. In section 4 we present our results. First we 
discuss the influence of the presence of fermions on the sphaleron. Then we exhibit 
the change in energy and shape of the barrier and present the fermion eigenvalue along 
the minimal energy path. We end this section with a discussion of the non-topological 
soliton, which, for heavy fermions, is reached at the end of the path instead of the 
vacuum. We present our conclusions in section 5. 



2 Weinberg-Salam lagrangian 

We consider the bosonic sector of the Weinberg-Salam theory in the limit of vanishing 
mixing angle. In this limit the U(l) field decouples and can consistently be set to zero 

C h = -\F; v F^ a + (D^(D^) - A($t$ _ 1^)2 (1) 
with the SU(2) L field strength tensor 

f% = W - + gt ab XK » ( 2 ) 

and the covariant derivative for the Higgs field 

D^={d ll -\igi»V*)*. (3) 

The SU(2)l gauge symmetry is spontaneously broken due to the non- vanishing vacuum 
expectation value v of the Higgs field 

<•»-;&(!)• (4) 

leading to the boson masses 

M w = M z = -gv , M H = vV2X . (5) 
2 
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We employ the values M\y = 80 GeV, g = 0.65. 

For vanishing mixing angle, considering only fermion doublets degenerate in mass, 
the fermion lagrangian reads 

- f (q) qL($u R + $rf R ) - f®(d R & + u R &)q L , (6) 

where q^ denotes the lefthanded doublet (ml, o?l), while q R abbreviates the righthanded 
singlets (u R , d R ), with covariant derivative 

D^ = (d,- l -igr a V;)q h , (7) 

and with <f> = ir 2 $*. The fermion mass is given by 

M F = . (8) 

Due to the U(l) anomaly baryon number and lepton number are not conserved 

d»j*' L = -hd»K» , (9) 

where 

= 167r2 -^"v- - ■ 3 - 



K » = T7r^W Tr (^ P V CT + \igV v V p V°) (10) 



{Fvp — ^-I^Flp, V<t — 1/2t 1 V^) is the Chern-Simons current and / g is the number 
of generations. In the unitary gauge the topological baryon number Q B , carried by a 
configuration, is determined by its Chern-Simons number N cs , 

N cs = J d 3 rK° . (11) 

For the vacua the Chern-Simons number is identical to the integer winding number, 
while the sphaleron at the top of the barrier carries half integer Chern-Simons number 
[5]. 



3 Equations of Motion 

The sphaleron barrier can be obtained by constructing a family of field configurations 
for the fermionic, gauge and Higgs fields, which interpolates smoothly from one vacuum 
sector to another, minimizing the energy along this path as a function of the Chern- 
Simons number [14]. Employing a spherically symmetric ansatz for the fields, we now 
derive the equations of motion, the energy functional and the Chern-Simons number. 
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3.1 Ansatz 

In the limit of vanishing mixing angle the general static, spherically symmetric ansatz 
for the gauge and Higgs fields is given by [22] 

„ = \^iM emifj+ m (Sm _ f ^ )+ fjM hf (12) 

gr gr gr 

V* = , * (13) 



$ = -^(H(r)+iT-rK(r))^ , 



(14) 



and involves the five radial functions /s( r )) fc{ r ), H(r) and K(r). 

To retain spherical symmetry we consider only fermion doublets degenerate in mass. 
The corresponding spherically symmetric ansatz for the fermion eigenstates is the 
hedgehog ansatz, 

q L (f,t) = e- lujt Ml[G L (r)+ta-fF L (r)} Xii , (15) 

q R (f,t) = e-^Ml[G R (r) - id ■ rF R (r)] Xh , (16) 
where the normalized hedgehog spinor Xh satisfies the spin-isospin relation 

<?Xh + rxh = 0. (17) 

The ansatz is form-invariant under a spherically symmetric gauge transformation 
with the unitary matrix 

U(r) = exp(i—^r-r) . (18) 

Under such a U(l) transformation the functions transform as follows 

J a + ifs — > exp(i<d)(f A + if B ) , 

fc — fc + r&, 
H + iK — ► exp(A (H + iK) , 

F L + iG L — > exp(i ® ) (F L + iG L ) , 

F R + iG R — > F R + iG R . (19) 
3.2 Energy functional 

The ansatz eqs.(12)-(17) leads to the spherically symmetric energy functional 

E = E b + E i , (20) 
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with the bosonic part 
4irM w 



E h 



9' 



dx 



r ^(/i + /l-i) 2 + (A + ^) 2 + (^ 



/a/i 



Cn2 



X 



X 



+ (K 2 + H 2 )(l + fi + f 2 B + J -±) + 2f A (K 2 



H 2 ) - Af B HK 



- 2xf c (K'H - KH') + 2x 2 {H' 2 + K' 2 ) + ex 2 (if 2 + if 2 - l) 2 



where 



4A 
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2^M^ 



) • 



(21) 
(22) 



and the fermionic part 



£f = 47rMi 



JO 



dxx" 



2 2 
F R G R — G R F R H — F R G R + F h G L — G L F L H — F L G L 



- 2 



G L F L 



— (G\ 

x 



Ft) 



2x K 1 



Ft) 



±LJ ^2x K ^ L 
+ 2M F H(G R G L - F R F L ) - 2M F K{F R G h + F L G R ) 

We have introduced the dimensionless coordinate 

x = M w r 

and expressed the fermion mass in units of M\y 

M F = M F /M W . 

The energy functional is invariant under the U(l) gauge transformation eq.(18). 

The fermion functions entering the energy functional need to be normalized. When 
N fermions occupy the eigenstate eqs.(15)-(16) the normalization condition is 



(23) 
(24) 
(25) 



roc 

An dxx 2 (G 2 R + F 2 + G 2 L + F 2 ) = N . 
Jo 



(26) 



To construct the minimal energy path over the sphaleron barrier, we then consider 
the functional 

W = E-u,N+*^ZN ca , (27) 

g 2 

where uj is the fermion eigenvalue and £ is a dimensionless lagrange multiplier [14]. The 
Chern-Simons number of a given configuration is 



cs 



-f 

2ir Jo 



dx 



(f 2 A + f 2 B )(--0')-( J( 



Jc 



, (28) 



. -e , )-( v / (/i + /l) sin(0-e))' 

where we have defined the function 6{x) = arctan(/ B //^). The function 9(x) is an 
arbitrary radial function, associated with a general U(l) gauge transformation (18). 
From the expression (28) the Chern-Simons number is readily obtained in an arbitrary 
gauge, and in particlular in the gauge, where the identification with the topological 
baryon number can be made. 
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3.3 Gauge choice 

We perform the calculations in the radial gauge, where fc = 0. In this gauge the spatial 
part of the Chern-Simons current contributes to the topological baryon number. We 
therefore rotate to the unitary gauge, where only the Chern-Simons number determines 
the topological baryon number. The corresponding gauge transformation involves the 
function 0(x), which satisfies 6(0) = and O(oo) = 9(oo). We then find for the 
Chern-Simons number the expression 

Ncs = ± jH dx(f B f' A - f A f B ) + ^ . (29) 

3.4 Equations 

Variation of the functional eq. (27) leads to the following set of equations for the bosons 

f'i = ^(fl + fl - 1) + Ia(K 2 + H 2 ) + K 2 -H 2 + g 2 xF L G L - U'b , (30) 

x z 

fs = %/l + /!-!) + fs(K 2 + H 2 ) - 2HK + \g 2 x{Gl - F*) + U'a , (31) 

Ob Zi 

^f B + e(H 2 + K 2 -l)H 

x z 

(32) 

K" = - 2 -K' + £- 2 ((l + fA f + f B ))-§f B + e(H 2 + K 2 -l)K 

- ^(F R G L + F L G R ) , (33) 
and for the fermions 

CuG L -F[- -F L + — — Fl - — G L + M F (-HG R + KF R ) = , (34) 

XXX 

CuF L + G' L + — — -^-Gl + — F L + M F (HF R + KG R ) = , (35) 

X X 

- CoG R + F R + —F R + M F (HG L - KF L ) = , (36) 

Ob 

wFr + G' R + Mf(HF l + KG L ) = . (37) 



+ - „ F (G R G L — F R F L ) , 
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Here u> is the fermion eigenvalue u in units of Mw 

a =ik- (38) 

For £ = 0, finite energy solutions of these equations correspond to extrema of the 
Weinberg-Salam theory, such as the sphaleron or the non-topological soliton. The above 
set of equations can then also be obtained by substituting the ansatz eqs.(12)-(17) into 
the general equations of motion. 

3.5 Boundary conditions 

The boundary conditions are chosen such that the energy density and the energy both 
are finite. 

At the origin the boson functions satisfy the boundary conditions 

f A (0) - 1 = f B (0) = H'(0) = K(0) = , (39) 
while the fermion functions satisfy 

Fr(0) = F L (0) = , (40) 

and 

G R (0) = c R , G L (0) = c L , (41) 

where c R , Cl are unknown constants, subject to the normalization condition (26). 
At infinity the gauge and Higgs field functions lie on a circle 

f A {oo) + i/ B (oo) = exp(i0(oo)) , H(oo) + itf(oo) = exp(*^) , (42) 

and 6{oo) is an unknown function of £. Therefore we choose the boundary conditions 

AH = f' B (oo) = H'(oo) = K'(oo) = (43) 
for the boson functions. The fermion functions all vanish at infinity 

F R (oo) = F L (oo) = G R (oo) = G L (oo) = . (44) 

4 Minimal Energy Path 

We construct the minimal energy path from one vacuum sector with iVcs = to the 
neighbouring vacuum sector with N cs = 1 by solving the equations of motion (30)-(37), 
varying the Chern-Simons number along the path by means of the lagrange multiplier 
£. For convenience we start at the top of the barrier, at the sphaleron solution, where 
£ = 0, and then vary the lagrange multiplier £ to positive and to negative values to 
follow the path down to the vacua on both sides of the barrier. 
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4.1 Sphaleron 



The sphaleron represents the configuration at the top of the barrier between topo- 
logically distinct vacua. It is a solution of the classical equations of motion of the 
Weinberg-Salam theory with one unstable mode. Previously the sphaleron solution 
was obtained only in the bosonic sector of the Weinberg-Salam theory, neglecting the 
fermions. Here we investigate the influence which the presence of fermions has on the 
sphaleron solution. 

Without fermions, the sphaleron (for 9 W = 0) has spherical symmetry and parity 
reflection symmetry. (We do not consider here the bisphalerons [23-24], where parity 
reflection symmetry is broken, and which arise at large values of the Higgs mass.) 
Parity reflection symmetry simplifies the ansatz for the bosonic fields of the sphaleron 
eqs.(12)-(14). Two of the three gauge field functions, namely /# and fc, and one of the 
Higgs field functions, H, vanish in the commonly used ansatz for the sphaleron. The 
sphaleron then has a Chern-Simons number of precisely 1/2 [5]. 

Fermions have a zero mode in the background field of the sphaleron. For zero mass 
fermions this zero mode is known analytically in terms of the sphaleron functions [11- 
13]. Recently this zero mode has been constructed explicitly for the case, where the 
fermions of a weak doublet are degenerate in mass [16]. In this case two of the four 
fermion functions, namely F& and F^, vanish. 

In a self-consistent treatment the presence of fermions affects the bosonic fields 
through the coupled equations of motion. The equation for the parity violating gauge 
field function fg now contains a source term from the fermions, which does not van- 
ish, even if F^ were zero. Similarly the equation for the previously vanishing Higgs 
field function H now contains a non-vanishing source term from the fermions (except 
for zero mass fermions, where Gr = 0). Thus the presence of fermions violates the 
parity reflection symmetry of the sphaleron. The previously vanishing parity violating 
functions now have source terms from the fermions and can therefore no longer vanish. 

As a consequence the fermions couple to a less symmetric bosonic configuration 
and therefore also the previously vanishing fermion functions F^ and F^ now have non- 
vanishing source terms, even if the fermion eigenvalue would still be zero. The self- 
consistent sphaleron solution thus involves eight radial functions, four for the bosons 
and four for the fermions. With the parity reflection symmetry of the sphaleron lost, 
the barrier no longer needs to be symmetric with respect to the Chern-Simons number 
iVcs — 1/2. Therefore the Chern-Simons number of the sphaleron at the top of the 
barrier may deviate from 1/2, and the eigenvalue of the fermions at the sphaleron 
solution may deviate from zero. 

In Fig. 1 we show the boson functions for the sphaleron for the Higgs mass Mg = 
Mw and for the fermion masses Mp = 130 GeV and Mp = 1.3 TeV for one bound 
fermion. We see, that even at a fermion mass of Mp = 130 GeV, corresponding to 
a mass on the order of the expected top quark mass, the functions fg and H are 
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still almost vanishing, and the deviation of the functions /a and K from those of 
the purely bosonic sphaleron is with only about 0.1% even an order of magnitude 
smaller. In contrast, for heavy fermions with masses in the TeV region f B and H differ 
considerably from zero, while the deviation of the functions j a an d K from those of 
the purely bosonic sphaleron is still small. 

In Fig. 2 we show the corresponding fermion functions. At the fermion mass of 
Mp = 130 GeV the functions and are still almost vanishing in the self-consistent 
calculation, they are smaller than 10~ 3 , while the functions Gr and Gl are almost 
identical to those of the background field calculation. For heavy fermions with masses 
in the TeV region the functions F-& and F^ differ markedly from zero. 

In Fig. 3 we show the energy of the sphaleron for three values of the Higgs mass, 
Mh = 50, 80, 100 GeV, as a function of the fermion mass, for one and three bound 
fermions, respectively. For massless fermions the energy of the self-consistent sphaleron 
is almost unaltered compared to the purely bosonic sphaleron. For Mh = M\y it 
increases by 0.0003 % (0.0034 %) for one (three) bound fermion(s). With increasing 
fermion mass the sphaleron energy first increases slightly until about Mp = 120 GeV 
and then decreases for larger values of the fermion mass. For a fermion mass on the 
order of the top quark mass, M F = 130 GeV, the sphaleron energy decreases by only 
0.003 % (0.028) with respect to zero mass fermions for M H = M w . While for fermion 
masses approaching the TeV region, the sphaleron energy decreases more strongly, 
e. g. for a fermion mass of 1.3 TeV the energy changes by 1.8 % (12.5 %) for one 
(three) bound fermion(s) and Mh = Mw- 

Fig. 4 shows the correponding energy in the boson fields for the sphaleron as a 
function of the fermion mass. Here we see an increase in energy for large fermion 
masses. The decrease of the total energy for large values of the fermions mass must 
thus be due to the fermion eigenvalue, which must decrease considerably for large values 
of the fermion mass. The fermion eigenvalue is shown in Fig. 5, where we observe indeed 
a strong decrease for heavy fermions. Thus the fermion eigenvalue of the self-consistent 
sphaleron deviates strongly from zero for heavy fermions. 

For light and intermediate fermions the eigenvalue is shown in Fig. 6. Here we 
observe small positive values of the eigenvalue. At a fermion mass of zero the eigenvalue 
starts at 0.07 GeV (0.2 GeV) for one (three) bound fermion(s), rises to a maximum of 
0.6 GeV (1.8 GeV) for a fermion mass on the order of the top quark mass, and then 
declines, passing zero for a fermion mass of about 0.2 TeV in both cases for M H = M w . 

In Fig. 7 the Chern-Simons number of the sphaleron solution is shown for the 
same set of Higgs masses and for one (three) bound fermion(s) as a function of the 
fermion mass. For zero mass fermions the Chern-Simons number approaches 0.499 
(0.497) for one (three) bound fermion(s). For fermions with mass M F = 130 GeV the 
corresponding Chern-Simons number is 0.495 (0.485), while for heavy fermions with 
mass M F = 1.3 TeV the Chern-Simons number changes by 10 % (31 %) for M H = M w . 
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Thus the Chern-Simons number differs appreciably from the value Nqs = 1/2 only for 
heavy fermions. 

4.2 Barrier 

When the barrier is calculated self-consistently in the presence of fermions the barrier 
is no longer symmetric with respect to Nqs — 1/2. The minimal energy path from the 
sphaleron down to the vacuum sector with Ncs = is shorter and steeper than the 
path down to the other side to the vacuum sector with iV cs = 1. 

The minimal energy path is shown in Fig. 8 for M# = My/ for one light fermion 
with Mp = 1/lOMw and for one (three) heavy fermion(s) with Mp = lOMw- For the 
light fermion the path agrees well with the purely bosonic path [14]. The deviation 
is very small at the top and increases slightly along the path towards both vacuum 
sectors. There the positive mass of the fermions is approached for Nqs — ► and the 
negative mass for iVcs — > 1, because there is now a free fermion present with positive 
or negative energy. 

For large values of the fermion mass, the barrier is altered considerably. With 
increasing fermion mass the top of the barrier moves to smaller values of the Chern- 
Simons number, making the barrier increasingly asymmetric. Beyond a fermion mass 
of about 0.7 GeV for one bound fermion the minimal energy path heading towards 
the Nqs = sector does not reach the bosonic vacuum configuration with one (three) 
free fermion(s) any more. Instead a new bound state is reached, a non-topological 
soliton [19-21], discussed below. The minimal energy path down to the other side of 
the sphaleron towards the N cs = 1 sector does reach a bosonic vacuum configuration 
with one (three) free fermion(s) with negative energy. 

In Fig. 9 we show the dependence of the Chern-Simons number on the lagrange 
multiplier £ along the minimal energy path for Mh = My/, for one light fermion with 
M F = l/10Mw and for one (three) heavy fermion(s) with M F = 10M W . For small 
fermion masses the relation is monotonic and almost linear. While Ncs varies from zero 
to one, the lagrange multiplier £ changes from -2 to +2. In contrast for heavy fermions, 
where non-topological soliton solutions exist, £ reaches a minimum close to -2, and 
then rapidly increases again, approaching zero at the non-topological soliton solution 
as shown in Fig. 10. (Note, that the Chern-Simons number of the non-topological 
soliton is nonzero.) 

In Fig. 11 we show the fermion eigenvalue oo along the minimal energy path in the 
self-consistent calculation for Mh = My/ for one light fermion with Mp = l/10My/ for 
one intermediate fermion with Mp = Mw, and for one (three) heavy fermion(s) with 
M F = 10M W . Comparing the eigenvalues obtained in the self-consistent calculation 
with those of the background field calculation [16] we see little change for light and 
intermediate fermions. The eigenvalue still looks (almost) antisymmetric with respect 
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to iVcs = 1/2. In contrast, for heavy fermions the self-consistent minimal energy path 
starts from the non-topological soliton and thus with an eigenvalue u < Mp. The 
eigenvalue decreases monotonically with increasing Chern-Simons number and reaches 
zero long before the sphaleron is reached (e. g. for Mp = 0.8 TeV the fermion mode has 
eigenvalue zero at Nqs = 0.4354, while the sphaleron is reached at Nqs = 0.4700 for 
Mn — Myy)- Then the eigenvalue gradually approaches the negative mass, uo — > —Mp. 

There exist critical values of the fermion mass where, at a given value of the Chern- 
Simons number Nqs, the bound state enters the positive or the negative continuum. For 
light and intermediate mass fermions the critical values of the self-consistent calculation 
almost coincide with those of the background field calculation [16]. For instance for 
Mp = 80 GeV the bound state enters the positive continuum at Nqs = 0.0773 and the 
negative continuum at Nqs = 0.9216 in the self-consistent calculation for one fermion, 
and at N cs = 0.0748 and N cs = 0.9252 in the background field calculation [16]. For 
heavy fermions the self-consistent solution reaches the non-topological soliton in the 
-^cs — sector, leading to a very different behaviour for Nqs — > 0. 

In Fig. 12 we show the value of the Chern-Simons number where, along the minimal 
energy path, the fermions have eigenvalue zero for Mh = My/ and one (three) bound 
fermion(s). For massless fermions the eigenvalue is zero precisely at a Chern-Simons 
number of 1/2. (The contributions to Nqs from the integral and from the angle 6(oo) in 
eq.(29) cancel to better than 10 -5 .) The Chern-Simons number of the zero mode then 
decreases with increasing fermion mass rapidly towards zero. The end configurations 
correspond to non-topological solitons with zero energy eigenstates, which are reached 
at Mp = 2.416 TeV (1.384 TeV) with N cs = 0.00038 (0.00038) for one (three) bound 
fermion(s). 

Let us now compare the self-consistent path with the path obtained in [19]. For very 
heavy fermions the self-consistent minimal energy path no longer decreases monotoni- 
cally from the sphaleron towards the vacua. Instead a new critical behaviour appears 
for fermion masses around 4 TeV. The path down to the vacuum sector N C s = turns 
backward at a first critical point and forward again at a second critical point. This is 
illustrated in Fig. 13 for the fermion masses M F = 4 TeV, M F = 6 TeV, and M F = 10 
TeV for Mh = My/. In Fig. 14 the correponding dependence of the lagrange multiplier 
£ on the Chern-Simons number iVcs is shown. While for Mp = 4 TeV and Mp = 6 
TeV there is only one configuration with £ = 0, there are three such configurations 
for M F = 10 TeV. Since all configurations with £ = correspond to solutions of the 
classical field equations of the (simplified) Weinberg-Salam theory, there are additional 
unstable solutions for these high fermions masses, new sphalerons. 

In [19] a restricted variational ansatz is used for the boson functions, which approach 
the vacuum configuration on both sides of the barrier, while the fermion functions are 
obtained perturbatively (in f^/g), and the fermion energy approaches the positive 
and negative free mass. Thus this variational calculation does not approach the non- 
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topological soliton. Not surprisingly this variational ansatz leads to a very different 
dependence of the barrier on the fermion mass. When the fermion mass runs across 
the TeV region critical values of the fermion mass are encountered, where a new local 
energy minimum first appears and then disappears again along path [19]. The new 
local minimum is associated with a new kind of soliton in [19], different from the non- 
topological soliton [19-21], and existing only for very heavy fermions with 9.0 TeV 
< M F < 12 A TeV for M H = M w . 

There are several reasons to doubt the argument for the existence of a new soliton, 
given in [19]. The calulation of [19] is only variational and therefore gives only an upper 
bound on the energy. Our calculation shows, that no such new smooth minimum ap- 
pears along the self-consistent path, although a different critical behaviour does appear 
for fermion masses around 4 TeV. Furthermore, for fermion masses in the critical range 
9.0 TeV < M F < 12.4 TeV the non-topological soliton has a negative energy, indicating 
the breakdown of the valence fermion approximation [20,21]. 

4.3 Non-topological soliton 

We have seen, that for heavy fermions one encounters a non-topological soliton along 
the minimal energy path when approaching the vacuum sector with Nqs = 0. We now 
take a closer look at this solution of the Weinberg-Salam model, which represents a 
minimum and not a saddle point like the sphaleron. 

Since the gauge field has little influence on the non-topological soliton [19], we first 
consider a simpler model by setting the gauge field equal to zero. In this case the 
equations for the left-handed and for the right-handed fermions are equivalent and we 
are essentially left with a chiral linear a-model, scaled up in energy since the vacuum 
expectation value of the Higgs field is about 25 times bigger than /„-. 

In Fig. 15 we show the total energy of the cr-model non-topological soliton as a 
function of the fermion mass for one bound fermion and for several values of the Higgs 
mass ranging from 50 GeV to 5 TeV. We see, that for a given Higgs mass, there exists 
a critical value of the fermion mass above which the fermion can form a stable bound 
state, the non-topological soliton, which is lower in energy than the mass of one free 
fermion. These results agree with those of [19], where a direct minimization technique 
was used. 

For a given fermion mass the energy of the non-topological soliton increases mono- 
tonically with the Higgs mass. Thus the energy is clearly bounded from above by the 
energy obtained in the non-linear cr-model, where the Higgs field is confined to the 
chiral circle [19-21]. The energy of the solitons becomes negative for fermion masses 
in the range 3.6 < M F < 3.8 TeV for the Higgs masses considered. This behaviour is 
largely due to the valence fermion approximation, since taking into account the effects 
of the Dirac sea leads to a qualitatively different behaviour [20-21]. 
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We observe that the soliton branches approach the free fermion branch in two dis- 
tinct ways, depending on the value of the Higgs mass. For low Higgs masses the energy 
along the soliton branch is always lower than the corresponding free energy, and the 
soliton branch simply bifurcates from the free fermion branch at a critical value of the 
fermion mass. Above a critical Higgs mass of about M# = 0.7 TeV for one bound 
fermion, however, the behaviour is different. Here the soliton branch first crosses the 
free branch, then reaches a cusp at a critical value of the fermion mass, where a new 
(unstable) branch emerges. The new branch then approaches the free branch from 
above. This latter behaviour is well known for non-topological solitons [25]. 

The different critical behaviour of the soliton branches at various Higgs masses is 
also seen in Fig. 16. There we show the fermion eigenvalue as a function of the fermion 
mass for the same set of values of the Higgs mass. When for large values of the Higgs 
mass the energy branch approaches the cusp, the fermion eigenvalue exhibits an infinite 
slope. 

In Fig. 17 we show the radial functions for the non-topological soliton for Mn = My/ 
and Mp = 2.6 TeV for one bound fermion. We see, that the fermions are localized in 
a small region of space, while the Higgs field approaches its vacuum expectation value, 
where H(x) — > 1 and K{x) — > 0, much more slowly. 

Finally we consider the effect of the presence of the gauge field on the non-topological 
solitons in the Weinberg-Salam model. In Fig. 18 we show the gauge field functions of 
the non-topological soliton for the physical gauge coupling g = 0.65 for Mh = My/ and 
Mp = 1.3 TeV and Mp = 2.6 TeV for one bound fermion. They differ only little from 
zero. 

The change in energy of the non-topological soliton due to the presence of the 
gauge field is smaller than 0.8 % for Mp = 2.6 TeV and for M# = My/. The deviation 
decreases for smaller fermion masses. The non-topological solitons possess a small, but 
finite Chern-Simons number, N cs < 0.0004 for Mp < 2.6 TeV and M H = M w . 

5 Conclusions 

In the Weinberg-Salam theory topologically distinct vacuum sectors are seperated by a 
barrier, whose height is determined by the sphaleron energy. Here we have constructed 
the minimal energy path over the barrier from one vacuum sector to the neighbouring 
one self-consistently in the valence fermion approximation. To retain spherical sym- 
metry we have neglected the weak mixing angle and assumed that the fermions of a 
doublet are degenerate in mass. While the first approximation is very good [17,18], the 
latter is badly broken and needs to be improved. 

The presence of the fermions affects the energy and the shape of the barrier. The 
sphaleron solution on top of the barrier loses its parity reflection symmetry due to the 
coupling to the fermions. It carries no longer exactly Chern-Simons number Nqs = 1/2, 
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and the fermion eigenvalue at the top of the barrier is no longer precisely zero. The 
deviations of the self-consistent sphaleron from the purely bosonic sphaleron depend 
on the mass and on the number of bound fermions. The presence of massless fermions 
hardly affects the sphaleron, while fermions with masses on the order of the top quark 
mass have a small but noticeable effect. The Chern-Simons number is decreased by 
about one percent from Ncs = 1/2, and the fermion eigenvalue differs slightly from 
zero, being on the order of one percent of the fermion mass. Though the height of 
the barrier is still insensitive in this mass range. Due to these small changes the 
barrier is very slightly tilted to the left, becoming asymmetric with respect to the top. 
The fermion eigenvalue along the self-consistent path differs little from the one of the 
background field calculation [16] for light and intermediate mass (top quark) fermions. 
We conclude that in this mass range the sphaleron and the barrier are well approximated 
by neglecting the fermions, and that the fermion eigenvalues are obtained with rather 
good accuracy in the background field calculation, as compared to the self-consistent 
valence fermion approximation. Clearly other approximations, such as neglecting the 
effect of the Dirac sea, still need to be investigated. 

For heavy fermions with masses approaching the TeV region the barrier deforms 
stronger. The height of the barrier then decreases monotonically with increasing 
fermion mass, and decreases the stronger the more fermions occupy the valence level. 
Likewise the Chern-Simons number of the sphaleron and the fermion eigenvalue at the 
sphaleron decrease. 

Heavy fermions with masses greater than about 0.7 TeV can form non-topological 
solitons in the vacuum sector. For such heavy fermions the self-consistent minimal 
energy path over the barrier does not end in the vacuum configuration but instead 
reaches the non-topological soliton as its end configuration. This is in contrast to 
the variational approximate path considered by Petriashvili [19], which for any value 
of the fermion mass ends in a vacuum with a free fermion. Due to this choice of end 
configuration Petriashvili finds a new minimum along the path for very heavy fermions, 
and therefore conjectures the existence of a new type of soliton [19]. Our self-consistent 
calculations also show a new critical behaviour for very heavy fermions. But we do not 
observe a new smooth minimum along the path. Instead we see the occurrence of 
two bifurcations for very heavy fermions. The minimal energy path then encounters 
two critical points. It winds backward at the first critical point, and forward again at 
the second critical point. This critical behaviour allows for the existence of additional 
unstable solutions, new sphalerons. Our self-consistent calculation does not support 
the conjecture, that a new type of soliton is present for very heavy fermions. 
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6 Figure Captions 



Figure 1: The gauge field functions Ja and fs and Higgs field functions H and K of 
the sphaleron are shown with respect to the dimensionless variable x for Mh = Mw 
and for M F = 130 GeV (solid) and Mp = 1.3 TeV (dotted). 

Figure 2: The righthanded and lefthanded fermion functions G and F of the 
sphaleron are shown with respect to the dimensionless variable x for M H = M w and 
for Mp = 130 GeV (solid and dot-dashed) and Mp = 1.3 TeV (dotted and dashed). 

Figure 3: The total energy (in TeV) of the sphaleron is shown as a function of the 
fermion mass (in TeV) for three values of the Higgs mass, Mh = 50 GeV, Mh = 80 
GeV, and M H = 100 GeV, for one fermion (solid) and three fermions (dotted). 

Figure 4: The boson energy (in TeV) of the sphaleron is shown as a function of 
the fermion mass (in TeV) for three values of the Higgs mass, M H = 50 GeV, M H = 80 
GeV, and M H = 100 GeV, for one fermion (solid) and three fermions (dotted). 

Figure 5: The fermion eigenvalue (in TeV) of the sphaleron is shown as a function 
of the fermion mass (in TeV) for three values of the Higgs mass, Mh = 50 GeV (solid), 
M H = 80 GeV (dashed), and M H = 100 GeV (dot-dashed), for one fermion and three 
fermions. 

Figure 6: The fermion eigenvalue (in GeV) of the sphaleron is shown as a function 
of the fermion mass (in TeV) for light and intermediate mass fermions for three values 
of the Higgs mass, M H = 50 GeV (solid), M H = 80 GeV (dashed), and M H = 100 GeV 
(dot-dashed), for one fermion and three fermions. 

Figure 7 : The Chern-Simons number of the sphaleron is shown as a function of 
the fermion mass (in TeV) for three values of the Higgs mass, Mh = 50 GeV (solid), 
M H = 80 GeV (dashed), and M H = 100 GeV (dot-dashed), for one fermion and three 
fermions. 

Figure 8: The energy (in TeV) along the minimal energy path is shown as a 
function of the Chern-Simons number for M H = M w for one fermion with Mp = 
1/10M W and Mp = 10M W (solid) and for three fermions with Mp = 10M W (dotted). 

Figure 9: The lagrange parameter £ along the minimal energy path is shown as 
a function of the Chern-Simons number for Mh = Mw for one fermion with Mp = 
l/10Mw and Mp = 10M W (solid) and for three fermions with Mp = 10M W (dotted). 

Figure 10: The lagrange parameter £ along the minimal energy path is shown in 
the vicinity of the non-topological soliton as a function of the Chern-Simons number 
for Mh = Mw for Mp = 10M^ for one (solid) and three (dotted) fermions. 

Figure 11: The eigenvalue (in units of Mp) along the minimal energy path is 
shown as a function of the Chern-Simons number for Mh = Mw for one fermion with 
Mp = 1/lOMw, M F = M w and Mp = 10M W (solid) and for three fermions with 
Mp = 10M W (dotted). 
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Figure 12: The Chern-Simons number of the zero mode along the minimal energy 
path is shown as a function of the fermion mass (in TeV) for Mh = Mw for one fermion 
and three fermions. 

Figure 13: The energy (in TeV) along the minimal energy path is shown as a 
function of the Chern-Simons number for Mh = My/ for one heavy fermion with 
Mp = 4 TeV, Mp = 6 TeV, and M F = 10 TeV. (The dotted parts are extrapolated.) 

Figure 14: The lagrange multiplier £ along the minimal energy path is shown 
as a function of the Chern-Simons number for Mh = Mw for one heavy fermion with 
Mp = 4 TeV Mp = 6 TeV, and Mp = 10 TeV. (The dotted parts contain non-negligible 
numerical errors.) The curves cross zero at sphaleron solutions. 

Figure 15: The energy (in TeV) of the non-topological soliton for vanishing gauge 
coupling constant is shown as a function of the fermion mass (in TeV) for Mh = 80 
GeV (dot-dashed), M H = 0.5 TeV (dashed), and M H = 1 TeV (dotted) and M H = 5 
TeV (solid) for one bound fermion. 

Figure 16: The fermion eigenvalue (in TeV) of the non-topological soliton is shown 
as a function of the fermion mass (in TeV) for Mh = 80 GeV (dot-dashed), Mh = 0.5 
TeV (dashed), and Mh = 1 TeV (dotted) and Mh = 5 TeV (solid) for one bound 
fermion for vanishing gauge coupling constant. 

Figure 17: The fermion field functions G (solid) and F (dotted) and Higgs field 
functions H (dashed) and K (dot-dashed) of the non-topological soliton are shown 
with respect to the dimensionless variable x for Mh = Mw and for Mp = 2.6 TeV for 
vanishing gauge coupling constant. 

Figure 18: The gauge field functions 1 — Ja and fs of the non-topological soliton 
are shown with respect to the dimensionless variable x for M H = M w and for M F = 2.6 
TeV (solid) and M F = 1.3 TeV (dotted) for the physical value of the gauge coupling. 
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